Acoustoelectric current and pumping in a ballistic quantum point contact 
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The acoustoelectric current induced by a surface acoustic wave (SAW) in a ballistic quantum 
point contact is considered using a quantum approach. We find that the current is of the "pumping" 
type and is not related to drag, i.e. to the momentum transfer from the wave to the electron gas. 
At gate voltages corresponding to the plateaus of the quantized conductance the current is small. It 
is peaked at the conductance step voltages. The peak current oscillates and decays with increasing 
SAW wavenumber for short wavelengths. These results contradict previous calculations, based on 
the classical Boltzmann equation. 
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The interaction of surface acoustic waves (SAW) with 
electrons in a two dimensional electron gas (2DEG) has 
recently attracted much attention. In particular the 
acoustoelectric effect (d.c. current driven by the SAW) 
was investigated experimentally in a point contact (PC) 
defined in GaAs/AlGaAs heterostructure by a split gate 
. Most of the theoretical considerations of this ef- 
fect were classical, based on the Boltzmann equation for 
electrons in a ID channel, with the SAW considered as 
a classical force Jl],£| or as a flux of monochromatic sur- 
face phonons d|j . Such an approach is valid only when 
the channel length is much longer than the electron Fermi 
wavelength and when the electron diffraction at the chan- 
nel ends can be neglected. In this picture the acousto- 
electric current results from the drag of electrons by the 
SAW. Its value is determined by the competition between 
the momentum transfer from the SAW to the 2DEG and 
the momentum relaxation due to impurity scattering [Q^] 
or due to electron escape from the PC §-§, for a ballistic 
PC. 

A quantum approach was used in [Q, but only PC's 
of length short compared to the SAW wavelength were 
considered for the experimentally relevant low frequen- 
cies. We present here a quantum description of the prob- 
lem, based on a different formalism, which allows a more 
general consideration and leads to results qualitatively 
different from those given by the classical approach. In 
particular, we find that the drag mechanism is not valid 
for the quantum acoustoelectric current, and that the 
reflection of the elctrons within the PC is crucial for pro- 
ducing the SAW effect. Unfortunately the results of the 
experiments do not allow to reach a definite conclusion 
about the mechanism of the acoustoelectric current. 

Consider a nanostructure (NS) of arbitrary geometry 



(e.g. a PC) where the 2DEG is confined by a poten- 
tial U(r) and is attached to terminals a (with no voltage 
bias). The NS is exposed to a random a.c. potential 
8U(r,t), produced by a gate or by radiation, infrared or 
acoustic, and localized within the NS. This a.c. potential 
induces a current through the NS, the d.c. component of 
which is the acoustoelectric current or the photovoltaic 
current, depending on the nature of the radiation. Using 
time- dependent scattering states (see below for details), 
we find that the d.c. current J/3 entering terminal f3 is 
given by (e < 0, % = 1) 



J/3 — ^ J/3a, Jf3a — — 

a 

dE ( df(E) 
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dE 



( 1 1 dridr 2 P(rx,T 2 ) x 
g (E\ ri ,T 3 )g a (E\ r 1; r 2 ). (1) 



The properties of the a.c. potential are condensed in the 
pumping factor 



^(ri,r 2 ) 



dwu8U(ri)8U(r 2 ) , 



(2) 



where 8U(ri)8U(r 2 ) is the Fourier component of the 
random field correlator 8U (ri, tx)8U (r 2 , t 2 ) and the over- 
bar denotes statistical averaging. The properties of the 
NS are embodied in 



Xan 



(3) 



where Xan(E\r) is a scattering state excited by an incom- 
ing wave w~ n (E\r) (normalized to unit incoming flux) of 
an electron with energy E entering the NS from channel 
n of terminal a. f(E) is the Fermi distribution in the 
terminals. 
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Equation (|l|) is valid in the weak field adiabatic approx- 
imation, when the a.c. perturbation SU is smaller than 
the Fermi energy Ep and when the relevant frequencies 
of this perturbation are smaller than all scales which de- 
termine the energy dependence of the scattering states. 
The statistical averaging replaces the temporal averag- 
ing, which is unavoidable when measuring a d.c. current 
induced by an a.c. perturbation. Below we assume zero 
temperature, which reduces the energy integration in (jl|) 
to E = Ep. 

The a.c. potential created by a SAW propagating in 
the x direction is SU(r,t) — A(t)expi(qx — u>ot) + c.c, 
where the amplitude A(t) is a stationary, slowly varying, 
random function. This potential is screened by the elec- 
trons of the 2DEG @. In the wide part of the PC the 
screening strongly reduces the potential (by the factor 
qap, where as is the Bohr radius), while in the narrow 
part screening is weak. To account for the screening effect 
we multiply SU (r, t) by a screening factor S(x) — S(—x). 
For this screened a.c. potential the pumping factor is 
P(ri, r 2 ) = 2ia;o|v4| 2 sing(xi - x 2 )S(x 1 )S(x 2 ) ■ 

Let us first consider a PC attached to two ideal ID 
leads at x — > ±00, and assume that the PC is represented 
by a repulsive delta function potential, U(x) — V5(x). 
The scattering states excited from the left terminal a = I 
at x — —00 and the right terminal a = r at x = +00 arc 
Xa{E\x) — v E {exp(±ikpx) + rp exp(ikp\x\)], where 
± denotes a = I and a = r, respectively, and kp 
and vp are the electron momentum and velocity at en- 
ergy E. The transmission and reflection amplitudes are 
tp = 1 + rp = (1 +%V/ve)~ 1 ■ For the screening factor we 
choose S(x) = exp(—\x\/ L s ). Eq. ([[]) then yields that 
the partial currents Ji r = J r i = 0, while Ju = —J rr = J, 
where J is the d.c. current through the PC in the in- 
direction. For q <C kp we have 



J 



uj \A\ 2 qki 



'2tt 2E 2 f q 2 
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where the index F means E — Ep. This result shows 
that (i) the current is proportional to wo, & n d hence is 
of the pumping type; (ii) the current increases with q for 
small wavevectors, and decays for large ones; (iii) a finite 
reflection is necessary for producing the effect. 

Turning now to a more realistic description of the PC, 
we model it by the 2D saddle-point potential U(x,y) — 
(l/2md 2 )[—(x/L) 2 + (y/d) 2 ], where m is the electron 
mass, L is the length of the PC and d is its width. 
For L 3> d this potential corresponds to a waveguide 
in the rr-direction with parabolic walls (at \x\ < L) ad- 
justed to horns (at |x| > L) with opening angle d/L. 
These horns represent the left and right terminals at x = 
q=oo. The scattering states are given by M Xan(E\r) — 
^n(y)x ± ( £ n\x). Here $„ is a normalized harmonic oscil- 
lator wave function with energy E n = A(n + 1/2), where 
A = 1/md 2 , and n = 0,1,2,... labels the channels in 



both terminals. (There is no channel mixing in a saddle- 
point potential.) X^{ E n\ x ) is given by the complex Weber 
(parabolic cylinder) function E(a,x), as defined in pj 

X ± (e n \x) = -iV^(Ld/2)^H(e n )E(~e n , ±0. (5) 

Here £ = (2/Ld) 1 / 2 x and s n = (E - E n )/5 with 5 = 
1/mLd. In Eq.(g) t(e) = (1 + e" 2 ^)" 1 / 2 is the transmis- 
sion amplitude of the barrier created by the saddle-point 
potential. Again, ± denote a = r, I, respectively. 

The Landauer conductance (at zero temperature, in 
units of e 2 /h) of such a PC is |11| Q = X) n ^( £ «) 2 , where 
now e„ = (Ep — E n )/ 5. When L> d the conductance as 
a function of Ep is a step like function, with plateaus of 
width A and steps of width S. The steps occur at energies 
E n where Ep equals the bottom of the transverse quan- 
tization mode n; For E > E n this mode is propagating, 
while at E < E n it is evanescent. 

The current in the saddle-point PC, as obtained from 
Eq. ([j]), consists of a sum over the separate mode con- 
tributions 



J = J ^F(e„,p), p = q(Ld/2) 1 / 2 , 



(6) 



with the nominal value J = 2e(aj Q /2Tr)(\A\ 2 /S 2 ). The 
function F(e,p) [see Eq. ( |l4| ) below] is positive for p > 0, 
and is odd in p, i.e. the electron flux is along the direc- 
tion of the SAW propagation. We find that F(e,p) is 
exponentially small when |e| 3> 1, that is, modes whose 
energies are far from the threshold E n do not contribute 
to the current. This is expected for the evanescent modes; 
for the propagating ones it means that in a free channel 
with no reflection the acoustoelectric current is zero. The 
crucial role of reflection in producing the current can be 
seen also from Eq. (0). 

Near the threshold, for |e| < 1, where the current is 
not small, we have (for p > 0) 



F{e,p) = 27re-' r£ t(e) 3 c(e)erf 



'2a 



P<1, 



F(e,p)=47rf(e) 2 cos 2 (L.-I- 7e 



P 



(7) 



Here a = Ld/L 2 , j £ — 2e\np + argT(l/2 — ie) and 
c(e) ~ 1 is given by integral H [see Eq. ( |l5|) below] at 
p = 0, cr = 0. In this calculation the screening factor is 
chosen to be S(x) — exp(— x 2 /L 2 ). Typically L s = L, 
which results in a = d/L <C 1. 

It follows from Eqs. (§) and (@) that for long SAW 
waves 5<L _1 the current increases linearly with q. For 
Cg< (id) -1 / 2 the current is independent of q and 
for short waves (Ld)^ 1 ^ 2 q it exhibits damped oscilla- 
tions. It is interesting to note that the oscillations are not 
simple geometrical; the wavelength "resonates" not with 
the channel length L, but with the less obvious length 
(Ld) 1 / 2 . The numerical calculations of a single mode 
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contribution to the current are depicted in the figure, for 
L/d = 10 and e = 0, ±0.5. The intermediate region in 
which F is independent of q is not distinguished for the 
not very small d/L ratios. One can see from this figure 
that below the threshold (e < 0) the current is much 
weaker then above it (e > 0). 



F(e,p) 
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Fig.l Contribution of a single channel to the current 
(see Eq. (6)) as a function of the dimensionless SAW 
wavenumber p, for different values of Fermi energy 
close to the step of the conductance 

Our quantum theory predicts that the current is strong 
only when Ep is at the vicinity of the transverse quan- 
tization channel bottom. This finding is in agreement 
with that of the classical approach Jl|,§-§), of giant os- 
cillations in the acoustoelectric current. However, in the 
quantum theory the width of the peak, S, is determined 
by the diffraction at the opening angle of the channel. 
Indeed, the width-to-spacing ratio of the oscillations is 
5/A = d/L. On the other hand, the peak width in the 
classical theory is determined (at zero temperature) by 
the scattering or by the escape time in the case where the 
SAW is described as a force, or by the energy and mo- 
mentum conservation when the SAW is described as a 
phonon flux. Experiment Q shows the oscillations, how- 
ever unfortunately the peaks are not very pronounced, 
and this is why one cannot make statements about the 
nature of their width. 

Another important difference between the quantum 
and the classical approaches concerns the behavior of 
the current at short waves. According to the (classi- 
cal) drag picture, the current should increase with q, 
whereas the quantum theory yields exponential decreas- 
ing. (Note that in these considerations the proportion- 
ality of ujo to q is not relevant.) For a short PC our 
theory predicts that J/ Jo — qL. This linear dependence 
on q was also obtained in Ref. [see Eq. (19) there], 
however the coefficient was not specified. Estimating in 
that equation cLT°(/i) ~ 1/5 (where T°(/j,) is the trans- 
mission at Ep ) the current (in our notation) becomes 
J/ Jo ~ qL x g(L/vsAw)$, where g is a geometry factor 
independent onL |. Since the ratio (L/vsaw)S is also 



i-independent, the predictions of both theories agree re- 
garding the dependences on SAW frequency and contact 
length. 

We now outline the derivation of Eq. ([I]). This is 
accomplished using the concept of time- dependent scat- 
tering states Q. Let the NS under the a.c. field be de- 
scribed by the Hamiltonian H = / dr*+(r)iJ(r, i)*(r), 
where H(r,t) = H (r) + SU(r,t) and H (r) = 
(l/277i)(-iV) 2 + [/(r). Here *(r) is the electron field op- 
erator. The time-dependent scattering state Xan (E\r, t) 
is defined as the solution of the equation 

i(d/dt) Xan (E\r, t) = H{r, t) Xan (E\r, t), (8) 

which is excited by an incoming wave w~ n (E\r) exp(-iEt) 
in the presence of the a.c. potential. The state 
Xan{E\r,t) is labeled according to the energy of the in- 
coming wave, but it contains components with energies 
E' =/= E, since due to the time dependent perturbation 
5U(r,t) the transmission and the reflection of the in- 
coming wave are inelastic. For a weak time-dependent 
potential, SU <C Ep, this equation can be solved by 
iterations, 

Xan (E\r,t)= (9) 
e-^[ Xan (E\r)+x^(E\r,t)+ X ^(E\r,t) + ...]. 

The first term here is the (time-independent) scattering 
solution of Hq, and the subsequent terms contain only 
outgoing waves. The latter can be found in terms of the 
retarded Green's function of Hq, 

(id/dt - flo(r))G(r, r',t) = S(r - r')S(t). (10) 

The time-dependent field operator, required for the 
calculation of the current density operator j(r, t) = 
— (ie/2m)^(r, i^V^r, t)+h.c, can be written in terms 
of the scattering states, 

/dE 
— J2^n(E) Xa n(E\r,t). (11) 
an 

Here a (E) is an operator creating an incoming electron 
in channel n of lead a with energy E. The averages of 
these operators are determined by the temperatures and 
the chemical potentials of the terminals connected to the 
leads. For the scattering states which are normalized to 
unit incoming flux 

(E')) = 2nS(E-E')6 an.a'n' ME), (12) 

where f a {E) is the Fermi distribution in terminal a. 

Using the above results one performs the quantum and 
statistical averaging to obtain the current density (j(r)). 
Evaluating (j(r)) far away in terminal (3 and integrating 
over the cross section of this terminal gives Jp of Eq. ([!]). 
The asymptotic behavior of the current density is derived 
using the following relation |T3| ] for the Fourier transform 
of the Green function defined by Eq. (|To|) , 
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G(E\t,t%^ = -iJ2™j m (E\*)X0m(E\r'), (13) 

771 

where wt m (E\r) is an outgoing wave in channel m of 
terminal (3 (normalized to unit flux). 

The acoustoelectric current for the saddle-point con- 
fining potential, Eq. (||), is obtained using for the com- 
plex Weber functions the representation jnj E(a, x) = 
k~ 1 / 2 W{a,x) + ik 1 ^ 2 W(a,—x), where W(a, ±x) are the 
real Weber functions and k = (1 + e 2 ^ ) 1 / 2 — e" a . We 
find 



F(e, P ) = t(e) 3 G(e, P )H(e,p), 



(14) 



where 



/+oo 
d£exp(-< 2 /2)sinKs(£,£), 
-oo 

/+oo 
^exp(-^ 2 /2)cosp^%,0, (15) 
-oo 



with 

g(e,0 = W(- 



-s,-Z) 2 -W(-e,0 2 = -g(e,-Q, 
= W(-e,£)W{-e,-Z)=h(e,-Q. 



(16) 



The appearance of the transmission amplitude t{e) in 
Eq. ([li]) ensures the exponential smallness of the func- 
tion F(e,p) for evanescent modes. To show that it is also 
exponentially small for propagating modes belonging to 
e > 1 we use the Darwin representation of the Weber 
functions (£ > 0) and obtain 



$(e,0 = [v^cosh(s/2)] [e-^ + sinfl], 

h(e,£) = [2Vecosh(s/2)]~ 1 cos6', 
6» = e(s + sinhs), £ = 2 v / esinh(s/2). 



(17) 



The functions g and h then contain exponentially small 
or fast oscillating terms. The integrals G and H can then 
be calculated near the saddle-point s = itt and are found 
to be ~ e-* e . 

We now turn to the case |e| < 1. For large p, F(e,p) 
is determined by the singular points of g and h. These 
functions are regular on the real £ axis; at £ — > 00 they 
are given by 



g(e,0 = 2C 1 [e- 7Te + t(er 1 sm$}, 
•d = £ 2 / 2 + 2eln£ + argT(l/2 - is). 



(18) 



Thus, the main contribution to H and G comes from 
large £, as both g and h have a singular point £ = 00. 
Using Eq. ( fDf ) one can check that the saddle points of 
cos ■& cos p£ and sin sin p£ are £ = ±p. Calculating G 
and if near these points yields the second of Eqs. ([?]). 

For small p the behavior of G and H is different. In 
if one can put cosp£ = 1 and a = 0. In G, however, the 



limit p — ► should be taken with care: for er = 0, G has 
a singularity of the form sgn(p) coming from the non- 
oscillating term in g. The factor erf(|> / \/2a) in the first 
of Eqs. (Q) results from the smoothing of this singularity 
by the screening factor. 

Being proportional to the frequency c^o of the SAW, 
the acoustoelectric current is of the pumping type: in- 
dependent of the value of u>o, a given fraction of the 
electron charge is transferred through the PC during 
each period of the SAW. Therefore, the current can be 
compared with the pumping current produced by two 
gates with phase shifted a.c. potentials, 5U(r,t) = 
Ai (t)ui (r) + A2(t)u2(r). Let the gates be symmet- 
ric, Ui(r) = u(x),u,2(r) = u{—x) and take A\{t) — 
A(t) cos(uj t + ip(t) + 0), A 2 (t) = A(t) cos(cj * + ¥>(*)), 
where the amplitude A{t) and the phase <p{t) are slowly 
varying, random functions, but the phase shift (f> is fixed. 

One can compare now the pumping factors, Eq. (||), 
for the SAW and the two gates and see that they 
are equal if one replaces exp(iqx) by [exp(?7r/4)u(a;) + 
exp(— wr/4)ii(— x)]/\/2- It means that a propagating 
SAW is equivalent to pumping with a phase shift <f> = n/2 
between two symmetric gates. This is why the acousto- 
electric current, although being of pumping type, does 
not contain the factor sin^>, typical for pumping pi 15 1. 
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